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SCATTER CORRECTION METHODS AND 
APPARATUS 

BACKGROUND OF THE INVENTION 

[0001] This invention relates generally to computed tomography and 
more particularly to cone beam volumetric computed tomography. 

[0002] An interference caused by x-ray scatter can undermine CT 
image quality and the validity of CT-based diagnostic imaging. In a single-row fan 
beam (FB) CT, the x-ray scatter is substantially reduced or almost eliminated through 
(1) shutting the aperture of a x-ray source collimator to a slit as the severity of x-ray 
scatter is proportional to the area illuminated by an x-ray beam, and (2) utilization of a 
detector collimator, which is cylindrically symmetric and focus at the focal spot of the 
x-ray source, to shield scattered x-ray photons but pass primary x-ray ones as much as 
possible. In a contemporary multi-row (up to 16) FB CT, those two measures are still 
being employed successfully to combat x-ray scatter interference. Hereinafter, both 
the single-row FB CT and multi-row FB CT are called FB CT. However, with a 
significant increment in row number, e.g., above 64, of a x-ray detector in cone beam 
(CB) volumetric CT, the techniques mentioned above will not work as well as their 
counterparts in FB CT, because (a) the aperture of x-ray source collimator is typically 
opened relatively wide in CB volumetric CT, and (b) the grid ratio of cylindrically 
symmetric detector collimator in volumetric CB CT can no longer be as high as that of 
their counterparts in FB CT, since an x-ray beam is intrinsically in sphere symmetry, 
and the uniformity of dynamic quantum efficiency (DQE) of an x-ray detector will be 
compromised by a cylindrically symmetric detector collimator if the grid ratio is as 
high as those in FB CT. 

BRIEF DESCRIPTION OF THE INVENTION 

[0003] In one aspect a method is provided. The method includes 
scanning an object with a cone beam volumetric computed tomography (CBVCT) 
system with a beam pass array positioned between an x-ray source of the CBVCT 
system and the object to acquire scatter data, scanning the object with the CBVCT 
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system without the beam pass array positioned between the x-ray source and the 
object to acquire image data, and correcting the image data using the scatter data. 

[0004] In another aspect, a cone beam volumetric computed 
tomography system is provided. The cone beam volumetric computed tomography 
system includes an x-ray source, a detector positioned to receive x-rays emitted from 
the source, and a beam pass array removably positioned between the source and the 
detector. 

[0005] In a further aspect, a computer readable medium encoded with 
a program is provided. The computer readable medium is configured to instruct a 
computer to receive scatter data from a cone beam scan of an object with a beam pass 
array present, receive image data from a cone beam scan the object without the beam 
pass array present; and correct the image data using the scatter data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0006] Figure 1 is a pictorial view of a CT imaging system. 

[0007] Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1 . 

[0008] Figure 3 illustrates a Bean Pass Array (BP A). 

[0009] Figure 4 illustrates the BP A shown in Figure 3 installed in the 
CT system shown in Figure 1 . 

DETAILED DESCRIPTION OF THE INVENTION 

[0010] Apparatus and methods based up on x-ray scatter intensity 
sampling and image sequence processing are herein described which combat x-ray 
scatter interference in CB volumetric CT. The apparatus and methods are illustrated 
with reference to the figures wherein similar numbers indicate the same elements in 
all figures. Such figures are intended to be illustrative rather than limiting and are 
included herewith to facilitate explanation of an exemplary embodiment of the 
apparatus and methods of the invention. 
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[0011] In some known CT imaging system configurations, an x-ray 
source projects a fan- shaped beam which is collimated to lie within an X-Y plane of a 
Cartesian coordinate system and generally referred to as an "imaging plane". The x- 
ray beam passes through an object being imaged, such as a patient. The beam, after 
being attenuated by the object, impinges upon an array of radiation detectors. The 
intensity of the attenuated radiation beam received at the detector array is dependent 
upon the attenuation of an x-ray beam by the object. Each detector element of the 
array produces a separate electrical signal that is a measurement of the beam intensity 
at the detector location. The intensity measurements from all the detectors are 
acquired separately to produce a transmission profile. 

[0012] In third generation CT systems, the x-ray source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged such that the angle at which the x-ray beam intersects the object 
constantly changes. A group of x-ray attenuation measurements, i.e., projection data, 
from the detector array at one gantry angle is referred to as a "view". A "scan" of the 
object comprises a set of views made at different, gantry angles, or view angles, during 
one revolution of the x-ray source and detector. 

[0013] In an axial scan, the projection data is processed to 
reconstruct an image that corresponds to a two dimensional slice taken through the 
object. One method for reconstructing an image from a set of projection data is 
referred to in the art as the filtered back projection technique. This process converts 
the attenuation measurements from a scan into integers called "CT numbers" or 
"Hounsfield units", which are used to control the brightness of a corresponding pixel 
on a cathode ray tube display. 

[0014] To reduce the total scan time, a "helical" scan may be 
performed. To perform a "helical" scan, the patient is moved while the data for the 
prescribed number of slices is acquired. Such a system generates a single helix from a 
fan beam helical scan. The helix mapped out by the fan beam yields projection data 
from which images in each prescribed slice may be reconstructed. 
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[0015] Reconstruction algorithms for helical scanning typically use 
helical weighing algorithms that weight the collected data as a function of view angle 
and detector channel index. Specifically, prior to a filtered backprojection process, 
the data is weighted according to a helical weighing factor, which is a function of both 
the gantry angle and detector angle. The weighted data is then processed to generate 
CT numbers and to construct an image that corresponds to a two dimensional slice 
taken through the object. 

[0016] To further improve the performance of the CT system, multi- 
slice CT systems are built. In such systems, multiple projections are acquired 
simultaneously with multiple detector rows. For example, by extending the x-ray 
width in the z-axis forming a cone beam, cone beam helical scans are performed. 
Similar to the case of helical scan, weighting functions are applied to the projection 
data prior to the filtered backprojection process. 

[0017] As used herein, an element or step recited in the singular and 
preceded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references 
to "one embodiment" of the present invention are not intended to be interpreted as 
excluding the existence of additional embodiments that also incorporate the recited 
features. 

[0018] Also as used herein, the phrase "reconstructing an image" is 
not intended to exclude embodiments of the present invention in which data 
representing an image is generated but a viewable image is not. Therefore, as used 
herein the term "image" broadly refers to both viewable images and data representing 
a viewable image. However, many embodiments generate (or are configured to 
generate) at least one viewable image. 

[0019] Referring to Figures 1 and 2, a multi-slice scanning imaging 
system, for example, a Computed Tomography (CT) imaging system 10, is shown as 
including a gantry 12 representative of a "third generation" CT imaging system. 
Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a detector 
array 18 on the opposite side of gantry 12. Detector array 18 is formed by a plurality 
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of detector rows (not shown) including a plurality of detector elements 20 which 
together sense the projected x-rays that pass through an object, such as a medical 
patient 22 between array 18 and source 14. Each detector element 20 produces an 
electrical signal that represents the intensity of an impinging x-ray beam and hence 
can be used to estimate the attenuation of the beam as it passes through object or 
patient 22. During a scan to acquire x-ray projection data, gantry 12 and the 
components mounted therein rotate about a center of rotation 24. Figure 2 shows only 
a single row of detector elements 20 (i.e., a detector row). However, multi-slice 
detector array 18 includes a plurality of parallel detector rows of detector elements 20 
such that projection data corresponding to a plurality of quasi-parallel or parallel slices 
can be acquired simultaneously during a scan to enable cone beam helical scans and 
therefore, system 10 is sometimes herein referred to as a cone beam volumetric 
computed tomography (CBVCT) system. 

[0020] Rotation of gantry 12 and the operation of x-ray source 14 are 
governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of 
gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 
analog data from detector elements 20 and converts the data to digital signals for 
subsequent processing. An image reconstructor 34 receives sampled and digitized x- 
ray data from DAS 32 and performs high-speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a storage device 38. 

[0021] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
cathode ray tube display 42 allows the operator to observe the reconstructed image 
and other data from computer 36. The operator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, x-ray 
controller 28 and gantry motor controller 30. In addition, computer 36 operates a 
table motor controller 44 which controls a motorized table 46 to position patient 22 in 
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gantry 12. Particularly, table 46 moves portions of patient 22 through gantry opening 
48. 

[0022] In one embodiment, computer 36 includes a device 50, for 
example, a floppy disk drive, CD-ROM drive, DVD drive, magnetic optical disk 
(MOD) device, or any other digital device including a network connecting device such 
as an Ethernet device for reading instructions and/or data from a computer-readable 
medium 52, such as a floppy disk, a CD-ROM, a DVD or an other digital source such 
as a network or the Internet, as well as yet to be developed digital means. In another 
embodiment, computer 36 executes instructions stored in firmware (not shown). 
Computer 36 is programmed to perform functions described herein, and as used 
herein, the term computer is not limited to just those integrated circuits referred to in 
the art as computers, but broadly refers to computers, processors, microcontrollers, 
microcomputers, programmable logic controllers, application specific integrated 
circuits, and other programmable circuits, and these terms are used interchangeably 
herein. 

[0023] Is has been well recognized that an x-ray projection image 
I t (x,y) consists of a primary image I p (x,y) formed by x-ray photons emanating 

from an x-ray source and a secondary image I s (jc, y) by scattered x-ray photons, 

I t (x 9 y) = I p (x,y) + I 5 (x,y) (1) 

and I s (x,y) can be further considered as a spatial low-pass filtering of I p (x,y) , 

I s (x,y) = LPF[I p (x 9 y)] (2) 

where LPF represents the spatial low-pass filtering. This means that the spatial 
variation of I s (x,y) is much lower than that of I p (x,y) . Hence, it is possible to 

reconstruct I s (x 9 y) from its samples measured at a sparse lattice. 

[0024] Figure 3 illustrates a Beam Pass Array (BPA) 60 including a 
plate 62 containing a plurality of holes 64 of an adequate diameter to remove or 
reduce x-ray scatter interference. In an exemplary embodiment, plate 62 is fabricated 
from lead. In one embodiment, each hole 62 is circular and has a diameter (d) of 
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about 1.5 mm, and spaced (D) from each other about 5 mm. Alternatively, holes 62 
are a shape other than circular such as a square where d refers to the length of a side. 
In another embodiment, diameter (d) is between 1.2 mm and 1.8 mm and space (D) is 
between 4 mm and 6 mm. In yet another embodiment, d is between 1.0 mm and 2.0 
mm and D between 3mm and 7 mm. 

[0025] Figure 4 illustrates BPA 60 removably positioned between x- 
ray source 14 and detector 18 in system 10. System 10 also includes a source 
collimator 70 and a detector collimator 72. A distance between x-ray source 14 and 
BPA 60 is designated as D S b. A distance between BPA 60 and the central axis of 
gantry 12 is designated Dbi. A distance between the central axis and detector 18 is 
designated as D jd . 

[0026] In use, two projection images are acquired. I x (x 9 y) without 
BPA 60 placed between x-ray source 14 and object 22 to be imaged, and I 2 (x, y) with 
BPA 60 placed between source 14 and object 22, in which x-ray intensities at 
illuminated spots are exclusively formed by primary x-ray photons, because the 
diameter of holes is adequately small. Subtracting I 2 (x 9 y) from I x (x 9 y) 9 one obtains 

h{x 9 y) = I x {x 9 y)-I 2 {x 9 y) (3) 

Measuring I A (x 9 y) at illuminated spots, one obtains I s (s x9 s ) 9 a sampled version 
ofI s (x 9 y) 

I s (s x9 s y ) = SS[I,(x 9 y)] (4) 

where the SS operator represents spatial sampling. Moreover, through spatial 
interpolation, such as linear or cubic B-spline interpolation, one can obtain an 
estimation of I 5 (x 9 y) 

I 5 (x,y) = SI[I s (s x ,s y )] (5) 

where the operator SI represents the spatial interpolation. Eventually, one obtains an 
estimation of the primary image 

7(x 9 y) = I l (x 9 y)-I s (x 9 y) (6) 
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in which x-ray scatter interference has been reduced or eliminated. 

[0027] A projection image sequence P t (a n ;x,y) (n = \...N) is 
generated by consecutively acquiring x-ray projection images along an x-ray source 
trajectory, where a n represents the angular position of x-ray source. Similarly, 

P t (a n ;x,y) (w = l..JV) can be decomposed into a primary image sequence and a 
scatter image sequence, wherein 

P t (a n ;x 9 y) = P p (a n ;x 9 y) + P s (a n ;x y y) (n = l...N) (7) 

[0028] According to the rationale elucidated above, two image 
sequences are acquired at each x-ray source position to obtain the primary image 
sequence P p (a n ;x 9 y) (n = l...N): (1) P x (a n ;x,y) (n = l...N) without a BPA; (2) 

P 2 (a n , ;x,y) (ri = 1 ...N') with a BPA placed between the x-ray source and the object. 

Ideally, a n should be equal to a n . , but in reality a n * a n , , because the angular position 

of x-ray source 14 may not be the identical among different gantry rotations. In one 
embodiment, the angular misalignment is corrected because one can obtain another 

image sequence P x (a n ,\x 9 y) (n' = l...N') that is aligned with P 2 {a n ,\x,y) 

(n' = l...N') in angular x-ray source positions by angularly interpolating image 

sequence P x (a n \x y y) (n = 1...AT), wherein 

PMn^y) = AI[P x {a n ^y)} («' = 1..J\T) (8) 

where the operator AI represents an angular (linear or cubic B-spline) interpolation. 
[0029] Subsequently, 

^K^^) = ^K^^)-^ 2 (^^^) («' = l.Jn (9) 
P s (a tt >;s x9 s y ) = SS[P A (a n .;x 9 y)] (n' = 1-JT) (10) 

[0030] Again, through angular interpolation and spatial interpolation 

respectively, 

PAa n ;s x9 s y )^AI[P s (a n .;s x9 s y )} (n = L..N) (11) 
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P s (a„,x,y) = SI[P s (a n ;s x ,s y )] (n = \...N) (12) 

one eventually obtains 

P p (a n ;x,y) = P l (a n ;x,y)-P s (a n ;x,y) (n = \...N) (13) 

_ f 

[0031] Tomographic images in which x-ray scatter interference has 
been significantly reduced or almost eliminated are reconstructed from the image 

sequence P p (a n ;x,y) through various known CB reconstruction algorithms. 
Additionally, although the image sequence P x (a n \x,y) (w = l..JV) was angularly 

interpolated to obtain P x {a n . ; jc, y) (n' = l...N') 9 it is contemplated that P 2 (a n .\x,y) 
{ri = 1...AT) can be likewise angularly interpolated to angularly align the image 
sequences. 

[0032] Typically, N' « N due to an intensive angular correlation in 
the scatter image sequence. The diameter (d) of holes 64 are adjusted appropriately to 
achieve an accurate measurement, since a too big hole can be passed through by 
scattered x-ray photons, and a too small hole is subject to the interference caused by 
penumbra of primary photons. Meanwhile, the inter-distance of (D) of holes 64 are 
adjusted appropriately to achieve an accurate estimation of the image formed by 
scattered x-ray photons, because a too sparse array can degrade the estimation 
accuracy of the scattered image and a too dense array can result in unnecessary x-ray 
dose to a patient 

. [0033] Additionally, although described in a medical setting, it is 
contemplated that the benefits of the invention accrue to all CT systems including 
industrial CT systems such as, for example, but not limited to, a baggage scanning CT 
system typically used in a transportation center such as, for example, but not limited 
to, an airport or a rail station. 

[0034] The herein described methods and apparatus improves image 
quality in CB volumetric CT by reducing or eliminating interference caused by x-ray 
scatter, and improves image quality in CB volumetric CT by maintaining a DQE 
uniformity of x-ray detectors with large numbers of rows, e.g., above 64. 



